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I.  INTRODUCTION 


The  spectrum  of  the  source  must  be  considered  when  r ctempting  to 
correct  observations  of  a source  for  the  effects  of  absorption  by  the  atm'*  * 
sphere  between  the  source  and  the  sensor.  Computer  programs  capable 
of  calculating  transmittance  and  radiance  over  an  inhomogeneous  path  at 
a spectral  resolution  fine  enough  to  include  individual  line  shapes  have  been 
developed  to  quantitatively  study  these  source-dependent  transmittance 
effects.  Because  these  computer  programs  are  applicable  to  a wide  variety 
of  additional  problems  involving  radiation  transport,  their  basis  and  use 
will  be  described  in  some  detail.  This  will  be  followed  by  a presentation 
of  results  computed  with  the  programs  applied  to  the  specific  problem  of 
quantitatively  estimating  the  effect  that  the  source  spectrum  can  have  on 
atmospheric  transmittance  corrections  applicable  to  measurements  from  a 
sensor  with  spectral  bandwidth  greater  than  a single  spectral  line  width. 


The  basic  problem  to  which  these  programs  have  been  applied  in  the 


present  study  is  summarized  in  Fig.  1.  A target  with  radiance  L>t  is  viewed 


by  a sensor  through  an  atmospheric  path.  The  idealized  sensor  indicates 


accurately  the  apparent  radiance  arriving  at  its  entrance  aperture.  A 


common  objective  of  such  measurements  is  to  infer  information  about  the 
source  radiance  from  the  observed  apparent  radiance.  To  achieve  this 
objective,  corrections  for  the  effects  of  the  intervening  atmosphere  must 
usually  be  made.  The  magnitude  of  these  corrections  depends  on  atmo- 
spheric conditions  along  the  path  from  source  to  sensor  and  the  spectral 
interval  over  which  the  sensor  operates.  Since  most  real  sensors  do  not 
have  a spectral  resolution  adequate  to  separate  individual  lines  in  the 
apparent  radiance  spectrum,  these  atmospheric  corrections  will  also  depend 
on  the  source  spectrum  and  the  spectral  resolution  of  the  sensor.  In  the 
infrared,  these  source-dependent  factors  can  be  particularly  significant  in 
cases  where  the  source  emission  is  due  to  hot  carbon  dioxide  and  water 


Preceding  page  blank 


to  the  Observation  of  Missiles  through  the  Atmosphere 


I 
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vapor,  since  these  same  gases  are  also  responsible  for  the  absorption  of 
radiation  by  the  atmosphere.  Of  course,  radiation  can  occur  from  sources 
other  than  hot  gases  and  be  prevented  from  reaching  the  sensor  by  atmo- 
spheric processes  other  than  molecular  absorption.  Some  of  these  other 
processes  are  listed  in  Fig.  1,  but  they  are  not  considered  in  detail  in  the 
present  study  since  they  are  not  expected  to  lead  to  significant  source- 
dependent  atmospheric  transmittance  corrections.  The  objects  of  the  studies 
reported  here  are  to:  (a)  demonstrate  the  magnitude  of  the  dependence  of 
atmospheric  transmittance  corrections  on  the  source  spectrum,  (b)  provide 
a means  of  computing  these  effects,  (c)  test  and  demonstrate-  the  use  of 
computer  programs  for  the  computation  of  high  resolution  transmittance 
and  radiance  spectra  for  inhomogeneous  paths,  and  (d)  test  the  validity  of 
the  AFCRL  line  atlas  for  hot  gas  spectral  computations. 

The  source  dependence  of  atmospheric  transmittance  corrections  for 
finite  spectral  bandwidth  sensors  can  be  computed  from  quasi-monochromatic 
source  radiance  and  atmospheric  transmittance  spectra.  The  unambiguous 
demonstration  of  such  effects  has  been  an  objective  of  the  present  study. 
However,  such  quasi-monochromatic  computations  are  much  too  time  con- 
suming for  the  routine  reduction  of  data.  For  this  reason,the  development 
of  a band  model  approach  to  the  efficient  computation  of  source -dependent 
transmittance  spectra  has  been  a part  of  the  present  study  and  is  reported 
in  Ref.  1.  A second  use  of  the  quasi-monochromatic  computations  has 
been  to  examine  the  validity  of  this  band  model  approach  to  source-dependent 
transmittance  computations.  This  comparison  is  presented  in  Ref.  1. 

The  computer  programs  for  radiation  transport  developed  as  part  of 
this  study  are  presented  in  some  detail  in  Section  II.  These  programs  are 
applicable  to  a wide  range  of  problems.  Results  obtained  with  these 


1 


S.  J.  Young,  Band  Model  Calculations  of  Atmospheric  Transmittance  for 

Hot  Gas  Line  Emission  Sources.  TR-0075(5647)- 1 , The  Aerospace 

Corporation,  El  Segundo,  California  (in  preparation,  June  1974). 
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computational  tools  are  given  in  Section  III.  Section  IV  is  devoted  to  a 
comparison  of  comjjuted  values  with  experimental  observations  to  ascertain 
the  applicability'  of  the  input  spectroscopic  data  and  the  accuracy  of  the 
programs . 
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II.  QUASI- MONOCHROMATIC  TRANSMITTANCE/ 
RADIANCE  COMPUTATION  METHOD 


A.  INTRODUCTION  AND  RADIATIVE  TRANSPORT 

EQUATIONS 

The  problem  of  particular  interest  in  this  stucy,  namely,  the 
computation  of  the  source  spectral  dependence  of  eflective  atmospheric 
transmittance,  is  a particular  type  of  radiative  transfer  problem.  In  study- 
ing such  problems  on  a quasi-monochromatic  basis  (as  is  done  in  this  study), 
complexities  arise  as  a result  of  he  masses  of  data  required  rather  than 
from  the  involved  mathematical  equations  that  describe  the  transport  of 
radiation.  The  basic  radiation  transport  equation  relating  for  a single 
frequency  v the  radiance  L (v)  at  location  0 to  the  radiance  L ( v ) at  the 
sensor  depends  on  the  transmittance  and  radiance  of  incremental  slabs  of 
the  intervening  atmosphere.  Referring  to  Fig.  2, 


L 
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(v 


) = Lt(v)T0D(v)+x; 


D 


(v)  <r.(w) 


(1) 


where 


i=n 

T „ = n t. 

m’  n i=m 


k(v,  z)  dz 


(2) 


These  equations  can  also  be  written  in  integral  form,  but  the  discrete  sum 
and  product  form  is  directly  applicable  to  implementation  on  a digital  com- 
puter and  so  will  be  retained.  Along  an  atmospheric  path  for  which  infrared 
absorption  is  significant,  it  is  usually  possible  to  assume  that  collisions 
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Fig.  2.  Physical  Relationships  Required  to  Develop 
the  Radiative  Transfer  Equation 
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between  molecules  occur  frequently  enough  for  a temperature  to  be  defined, 
and  thus  the  gas  is  in  local  thermodynamic  equilibrium.  The  temperature 
then  determines  the  distribution  of  energy  among  the  quantum  states  of  the 
molecule.  In  this  case,  the  radiance  of  an  atmospheric  segment  is 

<r.(v)  = [l  - t.(w)]  1^(1/,  0)  (3 

where  L^v,  0)  is  the  radiance  at  frequency  v of  a blackbody  at  temperature 
6.  For  this  case,  the  radiation  transport  problem  requires  only  the  com- 
putation of  the  transmittance  t.(w)  for  the  segments  along  the  optical  path. 

When  Eqs.  (1)  and  (3)  are  combined, 

N 

La(v)  = Lt(v)  T0,D(v)  + £ Tz. , D(v)  [l  - Ti(l/)J  Lb(v’e)  (4 

i-1  1 

Equations  (2)  and  (4)  form  the  basis  for  computer  calculation  of  the  trans- 
port of  radiatior  through  inhomogeneous  atmospheres.  The  Aerospace 
Corporation  radiation  transport  computer  program,  INHOM,  discussed  in 
more  detail  in  paragraph  II.  B employs  these  equations  to  approximate 
inhomogeneous  paths  by  up  to  81  homogeneous  segments. 

The  transmittance  of  individual  segments  can  be  computed  from 
spectroscopic  information  about  the  strength  and  width  of  the  spectral  lines 
resulting  from  the  various  possible  transitions  between  the  energy  levels 
of  the  molecules  in  the  oj  ical  path.  In  the  infrared  region,  these  tran- 
sitions are  normally  between  various  vibration- rotation  states.  There  are 
usually  a large  number  of  possible  transitions.  Only  with  the  recent 
appearance  of  the  compilation  of  line  parameters  obtained  under  the  direc- 
tion of  R.  A.  McClatchey  of  the  Air  Force  Cambridge  Research  Laboratory 


(AFCRL)2  has  there  been  available  the  necessary  spectroscopic  information 
organized  in  a fashion  to  facilitate  the  computation  of  atmospheric  spectra. 

The  AFCRi  line  atlas,  described  in  detail  in  Ref.  2,  provides 
spectroscopic  data  for  more  than  110,000  infrared  transitions  of  the  7 gases: 
HzO,  COz,  03,  N,,0,  CO,  CH4,  and  02<  For  each  transition,  the  strength 
S?  and  line  width  cfi  at  1-atm  pressure  and  296  K are  provided  along  with 
the  lower  energy  level  E?  of  the  transition  and  the  quantum  numbers  specify- 
ing the  energy  levels  between  which  the  transition  takes  place.  These  are 
arranged  in  th’  order  of  ascending  transition  frequency  v.  for  convenience 
in  spectral  calculations.  This  organization  mixes  the  transitions  of  various 
molecules  so  that  the  molecule  must  also  be  indicated  for  each  line. 

An  efficient,  but  flexible,  computer  program  operating  on  a fast 
digital  computer  is  required  to  utilize  any  spectroscopic  line  atlas,  such  as 
the  AFCRL  compilation,  for  calculations  covering  an  appreciable  spectral 
interval.  Our  approach  to  such  a computer  program,  designed  for  use  on 
the  Control  Data  7600  high  speed  computer  available  at  The  Aerospace 
Corporation,  is  described  in  the  remainder  of  this  section. 

B.  PROGRAM  DESIGN 

Equation  (4)  has  been  written  separating  the  source  radiance  from  the 
atmospheric  effects.  It  is  useful  to  maintain  this  separation  of  the  radiative 
transport  equation  into  source  and  atmospheric  terms,  as  is  implied  by 
Eq.  (4),  in  the  actual  computation  of  spectra  since  it  is  often  desired  to 
study  the  spectrum  resulting  from  the  transmittance  of  radiance  from  the 
same  source  through  a variety  of  optical  paths  or,  conversely,  the  trans- 
mittance of  various  sources  through  the  same  optical  path.  Tnis  separation 

2R.  A.  McClatchey,  W.  S.  Benedict,  S.  A.  Clough,  D.  E.  Burch, 

R.  F.  Calfee,  K.  Fox,  L.  S.  Rothman,  and  J.  S.  Garing,  AFCRL 
Atmospheric  Absorption  Line  Parameters  Compilation,  Environmental 
Research"Paper  434,  AFCRL- TR-73-0096,  Air  Force  Cambridge  Research 
Laboratories,  Mass.  (26  January  197  3). 
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then  becomes  desirable  because  the  combination  of  source  and  transmittance 
spectra  can  be  done  much  more  rapidly  than  the  computation  of  the  spectra. 
The  combination  operations  are  performed  by  a modified  version  of  an 
existing  program  that  also  has  the  capability  to  produce  plots  and  convolve 
instrument  functions  with  the  combined  spectra.  Figure  3 is  a schematic 
diagram  of  the  procedure  employed  to  compute  the  apparent  radiance  of  a 
target  after  the  radiation  has  passed  through  the  atmosphere.  In  this  pro- 
cedure, the  INHOM  Program  computes  quasi-monochromatic  spectra  based 
on  path  parameters  input  by  punched  cards  to  the  program.  It  obtains 
spectroscopic  data  from  the  AFCRL  line  atlas,  and  ‘he  resulting  spectra  are 
then  combined  by  the  SPECPST  Program.  The  capabilities  of  SPECPST 
are  listed  in  Table  1 and  described  in  detail  in  Ref.  3*  and,  therefore,  are 
not  delineated  in  this  report. 

Table  i.  SPECPST  Program  Capabilities 

• Combine  spectra  to  create  composite  spectra. 

• Convolve  spectra  with  defined  sensor  functions. 

• Compute  average  spectral  values. 

• Print  spectral  values  in  convenient  format. 

• Plot  spectra  for  desired  regions. 

Although  different  versions  of  the  INHOM  Program  are  shown  in  Fig.  3 
for  atmospheric  and  source  computations,  the  programs  differ  only  in  the 
subroutines  that  calculate  the  absorptance  of  the  segments.  The  relatively 
high  temperatures  and  absorber  concentrations  in  typical  plumes  require 
that  a number  of  factors  that  are  near  unity  for  atmospheric  spectra  must 


This  describes  the  basic  SPECPST  Program  only  and  does  not  include  the 
modifications  required  to  study  INHOM  produced  spectra. 

C.  M.  Randall,  User's  Guide  to  Computer  Programs  for  Fourier 
Spectroscopy,  TR  ^073(9260-01  )-8,  The  Aerospace  Corporation, 

El  Segundo,  California  (30  November  1972). 
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Fig.  3. 


Schematic  Diagram  of  the  Procedure  Employed  to  Compute 
the  Effect  of  the  Atmosphere  on  the  Transport  of  Infrared 
Radiation  (The  INHOM  Program  computes  the  quasi- 
monochromatic  transmittance  and  radiance.  The  SPECPST 
Program  has  the  capability  to:  combine  two  or  more 
spectra,  convolve  spectra  with  instrumental  functions,  and 
display  the  results  in  convenient  printed  and/or  plotted  form.  ) 
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be  specifically  included  in  the  source  spectrum  calculation.  A more 
detailed  description  of  the  INHOM  transmittance /radiance  program  and  its 
use  is  given  in  the  following  paragraphs. 

A radiative  transfer  program  should  be  flexible  to  allow  its  applica- 
tion to  a wide  variety  of  problems.  At  the  same  time,  however,  it  should 
be  fast  to  allow  economical  computing.  Satisfaction  of  these  two  somewhat 
contradictory  goals  in  the  INHOM  Program  is  attempted  by  separating  the 
various  required  types  of  calculations  into  subprogram  modules.  The 
modules  in  which  calculations  are  performed  that  are  not  1 _peated  often, 
e.g.  , the  transformation  of  atmospheric  parameters  into  parameters 
characterizing  the  various  segments  along  the  optical  path,  have  been  made 
very  flexible  and  general  with  only  minor  thought  given  to  computation 
efficiency.  On  the  other  hand,  those  modules  used  frequently,  e.g.,  to 
compute  the  line  shape  and  absorption  coefficient,  have  been  optimized  for 
efficient  execution,  even  to  the  extent  of  assembly  language  coding.  Flexi- 
bility for  these  computations  is  achieved  by  substituting  a different  version 
of  the  entire  module. 

The  basic  operation  of  the  INHOM  Program  is  indicated  schematically 
in  Fig.  4,  and  details  of  the  program  control  required  to  utilize  the  various 
options  are  discussed  in  the  Appendix.  Control  cards  are  read  giving  the 
program  information  until  a SPECTRUM  control  card  is  encountered  in  the 
input  stream.  This  card  initiates  the  computation  of  a spectrum.  The 
spectrum  computed  depends  on  information  communicated  on  the  SPECTRUM 
card  and  on  other  control  cards  encountered  previously  in  the  input  stream. 
The  optical  path  must  be  specified  prior  to  the  SPECTRUM  step  and  this 
may  be  done  in  two  ways.  It  may  be  specified  directly  by  giving  the  column 
density  of  each  gas  in  each  segment  along  with  the  temperature,  pressure, 
and  path  length  characterizing  that  segment,  or  it  may  be  specified  indirectly 
by  describing  an  atmosphere  in  terms  of  altitude  profiles  of  temperature, 
pressure,  and  varia  le  gas  concentration.  The  program  then  converts 
these  profiles  into  segment  parameters  by  employing  the  geometric  path 


Schematic  Diagram  of  INHOM  Program  Operation 


parameters  designated  by  the  user.  This  conversion  neglects  refraction 
effects,  and  assumes  a round  earth  with  a radius  of  6378.  16  km  surrounded 
by  symmetric  shells  of  atmos phe re . The  path  is  specified  by  giving  the 
zenith  angle  of  the  path  at  the  beginning  altitude  and  the  final  altitude.  Paths 
with  zenith  angles  greater  than  90  deg  are  acceptable.  If  the  final  altitude 
is  less  than  the  beginning  altitude,  two  paths  are  possible,  and  the  user 
may  select  either  one.  All  paths  that  intersect  the  earth's  surface  are 
recognized,  and  the  computation  of  such  a spectrum  is  prevented  and  an 
appropriate  error  message  is  printed.  Paths  entirely  outside  the  atmo- 
sphere are  also  recognized,  not  computed,  and  an  error  message  is  printed. 

The  conversion  of  an  inhomogeneous  path  to  a collection  of  homoge- 
neous segments  may  be  done  with  a number  of  different  assumptions.  Only 
a limited  amount  of  study  has  been  given  to  the  optimum  way  in  which  this 
should  be  done,  but  results  obtained  to  date  are  reported  in  paragraph  II. D. 

After  the  optical  path  is  properly  specified,  computation  of  the  radia- 
tion transport  equation  begins.  This  computation  is  indicated  by  the  DO 
statements  in  Fig.  4.  Computation  of  the  attenuation  in  a single  segment 
for  a single  frequency  is  done  in  a subroutine  called  SHAP  that  can  easily 
be  changed  to  accommodate  different  line- shape  functions. 

The  attenuation  k(v)  computed  by  the  SHAP  subroutine  for  a single 
frequency  v is  dependent  on  the  contribution  of  many  lines  in  the  neighborhood 

of  V. 


(5) 


where  IL  and  Iy  are  indices  bracketing  all  the  lines  assumed  to  be  close 
enough  to  v to  influence  k(v);  S.  is  the  strength  of  the  itn  line;  a.  is  the  width 


! 


of  the  ith  line;  f[(v  - v.),  a-.]  is  a line -shape  function;  a. 'id  W\  is  the  column 
density  of  the  gas  for  which  S.  and  <*.  characterize  a line.  The  factors  S. 
and  o.  are  dependent  on  the  pressure,  temperature,  and  species  concentra- 
tion in  the  particular  segment.  All  studies  to  date  have  been  carried  out 
assuming  a Lorentz  pressure-broadened  line-shape  function. 


fl(v  - V,),  a.N  ^ [(.  - v.)2  + »?] 


The  line  width  and  strength  parameters  depend  on  the  molecule  and  the 
conditions  in  the  segment  of  the  optical  path  being  considered.  The  AFCRL 
line  atlas  tabulates  the  line  strengths  and  widths  for  Pq  = 1 atm  and 
0q  = 296  K temperature.  Quantities  drawn  from  the  atlas  are  indicated  by 
a superscript  0.  In  most  of  the  calculations  reported  here,  line  strengths 
and  widths  appropriate  to  the  path  segments  have  been  obtained  by  the 
following  relations: 


'HE? 


s.  g(v,  0)  exp 


1 (®0  ' 


Qv  (e,m.)  Qr  (0,1m) 


Oj  = or?  (P/Po)(eo/0)1/2 


where  P and  0 are  the  average  pressure  and  temperature  in  a particular 
segment  of  the  optical  path,  and  h and  k are  the  Planck  constant  and  the 
Boltzmann  constant,  respectively.  The  correction  for  stimulated  emission 


g(v,0)  = 1 - exp(-hv/k0) 
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For  the  spectral  regions  presently  of  interest,  this  correction  differs 
significantly  from  unity  only  for  the  high  temperatures  involved  in  source 
radiance  computations.  The  term  is  the  energy  of  the  lower  energy  level 
of  the  transition,  and  (0,m.)  and  Qr  (9,m.)  are  respectively  the  ratio  of 
the  vibrational  and  rotational  partition  function  at  temperature  0 for  the 
molecule  rm  to  the  value  of  the  partition  function  at  0q  = 296  K.  For 
atmospheric  conditions,  differs  only  slightly  from  unity;  thus,  in  many  of 
the  calculation.-;  for  atmospheric  paths  reported  here  it  has  been  assumed  to 
be  unity,  ^e  rotational  partit'on  function  ratio  (0,rm)  is  approximated  by 

j_ 

Qr(0,m)  = (6/0 Q)  m (9) 

where  the  value  of  the  exponent  deprnds  on  the  molecule  being  considered. 

The  values  used  are  listed  in  Table  2 (Ref.  2).  For  source  conditions,  tho 

vibrational  partition  function  cannot  be  ignored  and  has  been  calculated 

4 

in  the  harmonic  approximation. 


-d. 

Qv(0,m)  = 77.  ll  - expf-hEj/kO)]  J (10) 

where  eY  is  the  energy  of  the  jth  fundamental  vibrational  level  and  cL  is  the 
degeneracy  associated  with  that  level.  The  coefficients  employed  are  listed 
in  Table  2.2’  4 

The  theory  of  line  broadening  is  complex  and  not  totally  adequate. 

For  atmospheric  conditions,  the  simple  approximation  of  Eq.  (7)  has  been 
employed.  However,  for  source  calculations,  the  concentration  of  certain 


4 

G.  Herzberg,  Molecular  Spectra  and  Molecular  Structure,  Van  Nostrand 
Re  inhold  Company,  New  York,  New  York  (1945),  p.  501 
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specie  may  be  high  enough  for  self-broadening  to  be  significant.  For 
these  cases,  we  have  used  the  results  of  Ludwig  et  al.  5 to  compute  the 
broadening  coefficient  to  be  applied  to  the  line  width  parameter  «?  from 
the  AFCRL  atlas  for  the  ith  line  of  species  m.  1,m 


l,  m 


°i,m  {L  y3m  Pj(0o/e)1/2  + Pm(e0/e)J/v^ 


(ii) 


where  is  the  broadening  coefficient  for  gas  j on  gas  m;  p.  is  the  partial 
pressure  of  gas  j,  is  the  resonant  self-broadening  coefficient;  and  y°  is 
the  result  of  carrying  out  the  computation  within  the  braces  for  dry  airTt 
eQ  = 296  K and  a total  pressure  of  PQ  = 1 atm.  These  are  the  conditions 


for  which  the  AFCRL  atlas  line  widths  are  calculated.  For  computational 

OPn nnrirta . iL re*  • i • q 

All  of  the 


economy,  the  coefficients  actually  used  are  /vC 

m Ym  Ym 


relevant  coefficients  based  on  Ref.  S are  lisle™  in  TrlbuY 
C'  CONTINUUM  ABSORPTION/TRANSMITTANCE 

All  absorption  mechanisms  that  vary  only  slightly  on  a spectral  seal, 
of  a few  line  widths  have  been  ie„ored  in  the  present  calculations,  since  the 
objective  of  the  present  study  i_  to  explore  the  effect,  of  line  correlation  on 
apparent  atmospheric  transmittance.  However,  capabilities  have  been 
incorporated  into  the  INHOM  Program  to  include  such  phenomena  when 
desired  bv  generating  an  appropriate  subroutine  module. 

The  choice  of  an  appropriate  bound  to  place  on  the  region  over  which 
lines  influence  the  absorptance  at  frequency  v,  I and  I in  Eq.  (5),  can  be 
viewed  as  a decision  on  where  to  begin  to  include  line  wings  as  a continuum 
phenomenon  to  he  handled  by  semiempirical  means  rather  than  by  considera- 
tion of  the  explicit  shape  of  individual  lines.  These  line-shape  functions  are 
5 — 

C.  B.  Ludwig,  W.  Malkmus,  J.  E.  Reardon  and  T a t 'ru 
Handbook  of  Infrared  Radiation  from  Cnmh„ ^ nasa™I!"’i 
Publication  NASA  SP-3080  (1973)7 — ^AbA  Special 
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Table  3.  Selective  Gas  Line- Broadening  Coefficients  for  H-O  and  CO 


probably  not  correct  at  great  distances  from  the  line  center.  An  appropriate 
bound  may  be  chosen  by  assuming  that  a Lorentz  line  shape  is  typical  and 
noting  that,  under  atmospheric  conditions,  Un'  v \dths  are  0.  1 cm  * or  less. 
The  rate  of  change  of  absorptance  then  decreases,  to  less  than  1 percent/ 
cm  beyond  10  cm  * from  the  line  center.  With  such  rates  of  absorptance 
change,  line  correlation  effects  will  be  minimal;  thus,  we  have  adopted  a 
limit  of  10  cm  away  from  a frequency  v for  the  inclusion  of  line  effects. 

D.  INHOMOGENEOUS  PATH  APPROXIMATION 

The  use  of  a series  and  continued  product  formulation  fox*  the  radiation 
transport  equation  in  an  inhomogeneous  medium,  as  defined  in  Eqs.  (2)  and 
(4),  implies  physically  that  the  continuously  varying  path  is  being  approxi- 
mated by  a series  of  homogeneous  segment?.  The  way  in  which  the  path  is 
divided  into  homogeneous  segments  and  the  number  of  segments  can  affect 
the  accuracy  of  the  computed  results.  As  a part  of  the  present  study,  a 
limited  variety  of  possible  path  division  algorithms  have  been  tried. 

Division  of  the  inhomogeneous  path  into  equal  geometric  length  seg- 
ments and  into  segments  containing  equal  numbers  of  air  molecules  has  been 
tried.  For  each  division  scheme,  a slant  path  through  a mid-latitude 
slimmer  model  atmosphere  from  0 to  20  km  at  a zenith  angle  of  60  deg  has 
been  divided  into  from  1 to  40  segments.  Then,  for  each  division  scheme 
and  each  number  of  segments,  the  transmittance  was  computed  in  t"  e 
15  cm  * band  from  3400  to  3415  cm  * with  a resolution  of  0.02  cm"  * . This 
spectral  region  was  chosen  because  it  contains  a range  of  transmittances 

_4 

from  0.  22  down  to  less  than  10  A spectrum  computed  for  40  segments 
is  shown  in  Fig.  5.  The  transmittance  spectra  resulting  from  various  runs 
are  compared  in  two  ways  in  Table  4.  First,  the  average  transmittance 
over  the  entire  15  cm  * band  for  the  various  number  of  segments  was  com- 
pared with  the  result  for  40  segments.  Second,  the  individual  transmittance 
spectra  were  divided  by  the  40- segment  transmittance  spectrum  for  all 
values  for  which  the  transmittance  was  greater  than  10”^.  In  the  typical 
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ratio  spectrum  shown  in  Fig.  6,  the  cross-hatched  region  indicates  where 
the  transmittance  was  less  than  1U  Peak  values  in  the  15  cm"1  band  for 
all  the  ratio  spectra  are  listed  in  the  columns  of  Table  4 headed  Peak  T / 
^40'  the  path  is  divided  into  40  segments,  the  division  method  makes 

little  difference,  but  for  a more  limited  number  of  segments,  the  equal 
number  of  air  molecules  yields  transmittances  that  are  closer  to  the  40- 
segment  values.  For  this  reason,  the  equal  absorber  concentration 
algorithm  has  been  adopted  for  the  calculations  reported  later  in  this  dis- 
cussion. When  calculations  over  inhomogeneous  paths  are  desired,  we 
have  chosen,  based  on  the  limicod  data  summarized  in  Table  4,  to  divide 
the  path  into  5 segments  This  choice  is  a compromise  between  the 
increased  accuracy  available  with  more  segments  and  the  increased  com- 
putation time  and  cost  required  to  carry  out  the  computations  over  a broad 
spectral  interval  with  an  increased  number  of  segments. 

E.  PROGRAM  EXECUTION  TIME 

The  time  t required  for  a spectrum  computation  is  directly  propor- 
tional to  the  product  of  the  number  of  segments  ng  into  which  the  optical  path 
is  divided  and  the  number  of  spectral  values  n^  desired  plus  a small  over- 
head time  tg  required  to  perform  those  calculations  needed  only  once  per 
spectrum. 


t = r n n + t_ 
s v 0 


(12) 


The  rate  constant  r depends  on  the  density  of  spectral  lines  in  the  spectral 
region  of  interest,  the  breadth  of  the  region  of  influence  included,  the  speed 
of  the  computer,  and  the  efficiency  of  the  program.  The  INHOM  Program 
described  in  this  report,  when  operated  on  the  CDC  7600  computer  at 
Aerospace,  achieved  a rate  of  r = 2.  2 msec  for  the  spectral  region  from 
3325  to  3975  cm"  and  an  influence  bound  of  10  cm"1.  Overhead  time  was 
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tQ  < 1 sec.  The  program  that  achieved  this  rate  included  all  the  pressure- 
broadening, stimulated-emission,  and  partition  function  calculations  described 
previously.  The  most  significant  single  improvement  in  speed  resulted  from 
coding  tiie  computer  operations  required  to  implement  Eq.  (5)  in  assembly 
language  with  an  ^.i-line  exponential  function  rather  than  in  FORTRAN.  This 
approach  resulted  in  a 25-percent  decrease  in  r. 
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III.  SOURCE  SPECTRAL  EFFECTS  ON 
APPARENT  ATMOSPHERIC  TRANSMITTANCE 

1 

I 


A.  TEST  CASES 

Based  on  initial  test  cases  described  in  the  following  paragraphs,  errors 
of  up  to  79  percent  are  possible  in  atmospheric  transmittance  corrections  if 
source  spectral  effects  are  not  considered.  The  initial  test  source  was  an 
isothermal  model  of  the  afterburning  plume  of  a typical  rocket  at  20  km. 
Characteristics  of  this  target  are  listed  in  Table  5.*  The  plume  is  assumed 
to  be  a right  circular  cylinder,  and  an  optical  path  along  a diameter  of  the 
cylinder  has  been  the  only  path  considered.  The  radiance  was  computed  at 
0.  01  cm  resolution  between  3325  and  3975  cm"1.  A short  segment  of  the 
target  spectrum  is  shewn  in  Fig.  7.  All  corrections  discussed  in  Section  II 
for  high  temperatures  were  included,  and  only  emissions  from  CC>2  and  H.,0 
were  considered.  No  lines  other  than  those  listed  in  the  AFCRL  atlas  were 
included.  This  approach  leads  to  an  error  of  unknown  magnitude  since  the 
AFCRL  atlas  was  developed  for  atmospheric  transmittance  work,  and  bands 
and  lines  arising  from  energy  levels  that  are  not  appreciably  populated  at 
atmospheric  temperatures  have  been  omitted  from  the  atlas.  Qualitative 
results  indica' o that  the  atlas  becomes  inadequate  for  temperatures  above 
about  1000  K.  Nevertheless,  we  have  used  the  atlas  to  gain  a feel  for  the 
magnitude  of  the  source  spectral  effects  on  atmospheric  transmittance, 
recognizing  tha.t  quantitative  results  may  not  be  correct. 

Studies  have  been  conducted  of  the  apparent  radiance  of  this  source 
after  passing  along  two  different  paths  through  the  tropical  atmosphere.  One 
is  a 100-km  horizontal  path  and  the  other  is  a path  to  space  at  a 75-deg  zenith 
angle  (see  Table  6).  A transmittance  spectrum  has  been  computed  for  each 
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Table  6.  Tropical  Model  Atmosphere^ 


Height 

(km) 

Pressure 

(mbar) 

Temp. 

(K) 

Density 

(g/m3) 

Water  Vapor 
(g/m3) 

Ozone 

(g/m3) 

0 

1.  013E+03 

300.  0 

1.  167E+03 

1.  9E+0« 

5.  6E-05 

1 

9.  040E+02 

294.  0 

1.  064E+03 

1.  3E+01 

5.  6E-05 

2 

8.  050E+02 

288.  0 

9.  689E+02 

? 3E+00 

5.4E-05 

3 

7.  150E+02 

284.  0 

8.  756E+02 

4.  7E+00 

5.  IE-05 

4 

6.  330E+02 

277.  0 

7.  951E+02 

2.  2E+00 

4.  7E-05 

5 

5.  590E+02 

270.  0 

7.  199E+02 

1.  5E+00 

4.  5E-05 

6 

4.  920E+02 

264.  0 

6.  501E+02 

8.  5E-01 

4.  3E-05 

7 

4.  320E+02 

257.  0 

5.  85  5E+02 

4.  7E-01 

4.  IE- 05 

8 

3.  780E+02 

250.  0 

5.  258E+02 

2.  5E-01 

3.  9E-05 

9 

3.  290E+02 

244.  0 

4.  708E+02 

.1.  2E-01 

3.  9E-05 

10 

2.  860E+02 

237.  0 

4.  202E+02 

5.  0E-02 

3.  9E-05 

11 

2.  470E+02 

230.  0 

3.  740E+02 

1.  7E-02 

4.  IE- 05 

12 

2.  130E+02 

224.  0 

3.  316E+02 

6.  OE-03 

4.  3E-05 

13 

1.  820E+02 

217.  0 

2.  929E+02 

1.  8E-03 

4.  5E-05 

14 

1.  560E+02 

210.  0 

2.  578E+02 

1.  0E-03 

4.  5E-05 

15 

1.  320E+02 

204.  0 

2.  260E+02 

7.  6E-04 

4.  7E-05 

16 

1.  110E+02 

197.  0 

1.  972E+02 

6.4E-04 

4.  7E-05 

17 

9.  370E+01 

i 9 5.  0 

1. 676E+02 

5.  6E-04 

6.  9E-05 

18 

7.  890E+01 

199.  0 

1.  382E+02 

5.  0E-04 

9.  OE-05 

19 

6.  660E+01 

203.  0 

1.  145E+02 

4.  9E-04 

1. 4E-04 

20 

5.  650E+01 

207.  0 

9.  515E+01 

4.  5E-04 

1.  9E-04 

21 

4.  800E+0! 

211.  0 

7.  938E+01 

5.  IE-04 

2.4E-04 

22 

4.  090E+01 

215.  0 

6.  645E+01 

5.  IE- 04 

2.  8E-04 

23 

3.  500E+01 

217.  0 

5. 618E+01 

5.  4E-04 

3.  2E-04 

24 

3.  OOOE+Ot 

219.  0 

4.  763E+01 

6.  OE-04 

3.  4E-04 

25 

2. 570E+0! 

221.  0 

4.  045E+01 

6.  7E-04 

3.  4E-04 

30 

1.  220E+01 

232.  0 

1.  831E+01 

3.  6E-04 

2.4E-04 

35 

6.  OOOE+OO 

243.  0 

8.  600E+00 

1.  IE- 04 

9.  2E-05 

40 

3.  050E+00 

254.  0 

4.  181E+00 

4.  3E-05 

4.  IE- 05 

45 

1.  590E+00 

265.  0 

2.  C97E+00 

1.  9E-05 

1.  3E-05 

50 

8.  540E-01 

270.  0 

1.  101E+00 

6.  3E-06 

4.  3E-06 

70 

5.  790E-02 

219.  0 

9.  210E-02 

1.  4E-07 

8.  6E-08 

100 

3.  000E-04 

210.  0 

5.  000E-04 

1.  0E-09 

4.  3E-11 

R.  A.  McClatchey,  R.  W.  Fenn,  J.  E.  A.  Selby,  F.  E.  Volz,  and  J.  S. 
Garing,  Optical  Properties  of  the  Atmosphere  (Third  Ed.  ),  Environmental 
Research  Paper  No.  411,  AFCRL-72-0497,  Air  Force  Cambridge  Research 
Laboratories,  Mass.  (24  August  1972). 
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of  these  paths  at  0.  01  cm"*  resolution  from  3325  to  3975  cm"*.  A sample 
segment  of  such  a spectrum  is  illustrated  in  Fig.  8.  In  addition  -.n  apparent 
radiance  spectrum  has  been  obtained  for  each  of  these  paths  by  rrultiplying 
the  transmittance  spectrum  for  the  path  by  the  source  radiance  function.  A 
sample  segment  of  such  an  apparent  radiance  spectrum  for  the  75-deg,  zenith 
angle  slant  path  is  shown  in  Fig.  9. 

B.  ATMOSPHERIC  TRANSMITTANCE  EFFECTIVE  FOR  LOW 
RESOLUTION  SENSOR 

A variety  of  atmospheric  transmittance  values  can  be  obtained  from  the 
spectra  presented  in  paragraph  III.  A when  a sensor  with  a spectral  bandwidth 
greater  than  a few  line  widths  is  considered.  The  differences  depend  on  the 
way  in  which  the  low  spectral  resolution  transmittance  is  defined. 

The  effective  average  transmittance  Tg  is  th**  plicable  value  if  the 
objective  of  sensor  measurements  is  the  inference  ol  the  source  radiance. 
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is  the  ratio  of  the  integrated  apparent  radiance  L (v)  to  the  integrated  target 


radiance  L.(\>).  Spectra  of  T for  the  two  optical  paths  considered  are  indicated 

t 6 , 


by  solid  lines  in  Figs.  10  and  11.  The  value  Av 


- Vj  = 20  cm 


-1 


was 


used  in  computing  these  spectra.  Three  values  of  T obtained  for  bandwidths 

- 1 e 

greater  than  20  cm  in  the  two  wings  and  the  center  of  the  absorption  region 


shown  in  Figs.  10  and  11  are  listed  in  Table  7.  The  entry  in  Table  7 titled 
"20  cm  * Intermediate  Average"  is  the  average  of  the  T curve  from  Fig.  10 


or  1 1 for  the  same  radiometer  bandwidth.  This  indicates  that  the  average 
of  low  resolution  spectra  does  not  produce  the  same  results  as  the  average 
obtained  from  high  resolution  spectra. 
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Transmittance  of  the  Atmosphere  along  a Path  from  20  km  to 
Space  at  a 75-Deg  Zenith  Angle  through  a Tropical  Atmosphere 


Fig.  9.  Apparent  Spectral  Radiance  of  a Typical 
through  the  Atmosphere  from  ZO  km  to  S 


FREQUENCY 


Table  7.  Radiometer  Transmittance 


; 


The  average  transmittance  T is  the  most  frequently  computed  or  ap- 
proximated transmittance. 


T =Sv/  2 T(v)dw  (14) 

V1 

The  average  transmittance  for  the  two  sample  paths  was  computed  from  the 
high  resolution  transmittance  spectra  by  averaging  over  Av  = 20  cm"*,  Re- 
sults are  indicated  by  dashed  curves  in  Figs.  10  and  11.  Average  transmit- 
tances  for  wider  bands  are  shown  in  Table  7.  For  sources  of  uniform 

spectral  radiance,  T should  be  the  same  as  T , but  for  line  sources  such  as 

e 

hot  gases,  they  are  clearly  not  the  same. 

The  factor  T is  the  quantity  approximated  by  most  atmospheric  trans- 
mittance band  model  procedures.  Estimates  of  T by  two  such  programs  for 
the  two  sample  paths  are  included  for  comparison  in  Table  7.  In  the  table, 

LOWTRAN  refers  to  the  program  developed  by  R.  A.  McClatchey  at  the  Air 

6 7 

Force  Cambridge  Research  Laboratories,  ’ and  ATLES  refers  to  a band 
model  developed  at  Aerospace  to  account  for  source  effects.  * Spectra  com- 
puted with  the  use  of  LOWTRAN  are  indicated  by  the  broken  curves  in  Figs. 
10  and  11.  The  comparisons  given  in  Table  7 indicate  that  the  band  models 
approximate  T reasonably  well.  The  discrepancy  between  the  T spectrum 
and  the  LOWTRAN  estimate  for  the  75-deg  slant  path  is  not  understood  r t 
present. 

An  attempt  to  include  source  effects  can  be  made  by  using  low  resolu- 
tion source  and  transmittance  functions  to  compute  low  resolution  effective 

transmittance  T T . 

eL 


J.  E.  A.  Selby  and  R.  A.  McClatchey,  Atmospheric  Transmittance  from 
0.  25  to  28.  5 M-m:  Computer  Code  LOWTRAN2,  Environmental  Research 
Paper  No.  <27,  AFCRL-72-0745,  Air  Force  Cambridge  Research  Labora- 
tories, Mass.  (29  December  1972). 


; 


; 

: 
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The  results  of  carrying  this  out  for  the  wings  and  center  of  the  2.  7-pm  absorp- 
tion are  given  in  Table  7.  The  accuracy  of  estimating  Tg  in  this  way  is  only 
a little  better  than  the  average  transmittance  T obtained  without  considering 
the  source  spectrum  in  any  way. 

The  techniques  of  band  modeling  can  be  applied  directly  to  the  definition 

of  the  effective  average  transmittance,  Eq.  (13),  to  obtain  Y These  pro- 

e o 

cedures  are  discussed  in  detail  in  Ref.  1,  as  are  spectral  comparisons  of  the 
effective  average  transmittances  derived  from  the  band  model  and  from  high 
resolution  calculations.  Included  in  Table  7 is  the  radiometer  transmittance 
predicted  by  this  special  band  model,  which  is  the  only  procedure  thr.t  explic- 
itly takes  into  account  line-correlation  relations  between  the  source  and 
atmosphere  in  estimating  T^.  It  is  also  the  only  procedure  that  does  not 
overestimate  the  transmittance.  In  fact,  if  any  conclusions  may  possibly  be 
drawn  about  the  band  model  from  this  limited  set  of  data,  it  would  appear 
that  the  model  overestimates  the  effects  of  line  correlation. 

Based  on  the  limited  test  cases  presented  in  this  discussion,  one  may 
conclude  that  all  average  transmittance  procedures,  which  do  not  account  for 
possible  line  correlations,  tend  to  overestimate  the  effective  transmittance 
of  the  atmosphere  for  a line  source.  In  some  spectral  regions  and  over  some 
atmospheric  paths,  this  can  lead  to  errors  of  up  to  79  percent  in  the  estimated 
effective  transmittance. 


i 


3 
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IV.  PROGRAM  AND  LINE  ATLAS  EVALUATION 


A.  COMPARISON  OF  COMPUTED  SPECTRA  AND  EXPERIMENTAL 

RESULTS 

A comparison  of  computed  and  experimentally  determined  spectra  has 
been  made  for  a limited  number  of  cases  to  demonstrate  the  overall  accuracy 
of  the  high  resolution  computation  procedures  outlined  when  they  are  applied 
with  the  AFCRL  line  atlas.  Only  H20  and  CO^  have  been  studied  as  they  are 
the  principal  atmospheric  absorbers  and  source  emitters  in  the  spectral 
region  of  interest  in  the  present  study.  Four  experimental  spectra  were 
chosen,  two  characteristic  of  atmospheric  paths  and  two  somewhat  charac- 
teristic of  missile  plumes.  The  conditions  for  each  case  are  listed  in  Table 
8,  and  the  spectra  are  illustrated  in  Figs.  12  through  16.  Experimental 
spectra  are  indicated  by  solid  lines  and  calculated  spectra  by  dashed  lines 
in  the  figures. 

B.  ATMOSPHERIC  COz  PATH 

Experimental  C02  spectra  characteristic  of  atmospheric  paths  have 
been  obtained  from  the  work  of  Burch,  and  his  Sample  10  has  been  chosen 
lor  detailed  comparison.  A low  resolution  spectrum  for  this  sample,  il- 
lustrated in  Fig.  12  by  a solid  line,  was  obtained  by  differentiating  the  inte- 
grated absorption  spectrum  tabulated  in  Ref.  8 and  converting  the  result  to 
a transmittance  spectrum.  The  computed  comparison  spectrum  (dashed  line 
in  Fig.  1L)  was  derived  by  convolving  a high  resolution  spectrum  obtained  by 
means  of  the  INHOM  Program  from  the  AFCRL  line  atlas  with  a 2.  5 cm’* 
full  width  at  half  maximum  (FWHM)  triangular  instrument  function.  Except 
in  the  region  from  about  3660  to  3725  cn.'*  agreement  of  the  calculated  and 


D.  E.  Burch,  D.  A.  Gryvnak,  and  R K.  Patty,  Absorption  by  CO?  Between 
3100  and  4100  cm~  (2.  44  - 3.  22  Mk  ns),  U-4132,  Aeronutronic  Division 
of  Philco-Ford,  Newport  Beach,  C.  1 , -nia  (30  April  1968). 
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Table  8.  Conditions  for  Experimental  Comparison  with  Computed  Spectra 


oratories  (August  1965)7 
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experimental  values  is  good.  In  the  region  where  agreement  is  less 
satisfactory,  the  apparent  disagreement  may  be  the  result  of  the  procedures 
employed  to  obtain  the  spectrum  plotted  as  the  experimental  spectrum.  How- 
ever, the  high  resolution  comparison  made  for  these  same  experimental 
conditions,  discussed  in  the  following  paragraphs,  indicates  there  is  a real 
difference  between  the  calculated  and  experimental  values  in  this  spectral 
region. 

Figure  13  is  a higher  resolution  spectrum  for  the  same  CC?  sample  of 

g ^ 

l»urch.  The  experimental  spectrum  was  obtained  in  this  case  by  photo- 
graphing the  spectrum  in  Ref.  8,  enlarging  it,  and  then  using  an  "OSCAR" 
digitizing  macht  « to  punch  on  cards  the  x and  y- coordinates  of  minima,  max- 
ima, and  inflection  points  that  characterize  the  spectrum.  These  cards  con- 
taining x and  y coordinates  then  became  input  data  to  a computer  program  to 
produce  the  linear  wave  number  plot  indicated  by  the  solid  curve  in  Fig.  13. 
When  the  initial  computer- generated  plot  of  the  experimental  data  was  com- 
pared with  the  plot  of  the  calculated  spectrum,  there  appeared  to  be  a con- 
stant shift  of  major  spectral  features,  though  the  major  features  seemed  to 
be  the  correct  distance  apart  spectrally.  Since  the  entire  experimental 
frequency  scale  depends  on  the  accuracy  of  the  location  of  two  tick  marks 
placed  on  the  original  spectrum  in  Ref.  8,  we  felt  justified  in  slipping  the 
experimental  spectra  in  frequency  space  until  the  major  features  agreed  in 
location  with  the  calculated  features.  This  required  a shift  of  less  than 
1 cm'*.  Nevertheless,  for  this  reason  absolute  frequency  comparisons  are 
not  possible  from  Fig.  13.  The  detailed  shape  of  individual  spectral  features 
may  also  be  distorted  slightly  in  Fig.  13  because  of  finite  sampling  of  the  .*pec 
trum  in  the  digitization  process.  The  computed  spectrum  indicated  by  the 
dashed  line  was  obtained  by  convolving  a triangular  instrument  function 
0.  8 cm’*  wide  (FWHM)  with  the  high  resolution  spectrum  obtained  for  the 
experimental  conditions  from  the  AFCRL  line  atlas  by  means  of  the  INHOM 
computer  program. 


j 
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Fig.  13.  Comparison  of  Computed  CO2  Transmittance  with  High  Resolution  Laboratory  Experimental 
Observations  Characteristic  of  Atmospheric  Paths  (The  conditions  are  listed  in  Table  8. 
The  solid  curve  was  obtained  from  the  spectrum  published  in  Ref.  8 by  means  discussed  in 
the  text.  The  dashed  curve  was  calculated  from  the  AFCRL  atlas  of  spectroscopic  para- 
meters by  means  of  the  computer  programs  discussed  in  this  report.  ) 


r 

! ! 

Based  on  the  CC>2  spectra  for  a single  sample  presented  here,  the 
AFCRL  line  atlas  and  INHOM  Program  appear  capable  of  predicting  with 
adequate  accuracy  the  absorption  of  CC>2  in  the  2.  7-pm  spectral  region  for 
typical  terrestrial  atmospheric  conditions.  Agreement  between  the  computed 
and  experimental  spectra  in  Fig.  13  is  good  except  in  the  region  from  3660 
to  3720  cm"*,  where  the  computed  transmittance  is  somewhat  higher  than 
the  experimental  transmittance,  and  below  about  3520  cm  *,  where  the 
positions  of  various  features  do  not  seem  to  agree  very  well.  The  difference 
between  3660  and  3720  cm"  1 is  probably  real,  while  the  low  wave  number 
difference  is  probably  uncorrected  scale  shift  in  our  reproduction  of  the 
original  data  and  is  probably  not  real. 

C.  ATMOSPHERIC  H20  PATH 

_ <j 

Figure  14  is  a high  resolution  spectrum  of  water  obtained  by  Burch 
and  identified  as  Sample  39.  The  experimental  transmittance  spectrum  was 
tabulated  in  Ref.  9,  and  thus  no  reproduction  scale  changes  or  interpolations 
have  been  required.  The  computed  spectrum  shown  by  the  dashed  curve  in 
Fig.  14  was  obtained  by  convolving  the  triangular  instrument  function  0.  5 
cm"1  wide  (FWHM)  with  the  high  resolution  spectrum  obtained  from  the 
AFCRL  line  atlas  by  means  of  the  INHOM  computer  program,  which  em- 
ployed the  experimental  conditions  listed  in  Table  8. 

Throughout  most  of  the  spectrum,  agreement  between  computed  and 
experimental  values  is  good.  However,  at  the  low  wave  number  end  of  the 
spectrum,  there  are  a number  of  spectral  features  in  the  computed  spectrum 
that  do  not  appear  in  the  experimental  spectrum.  These  have  b.  en  traced  to 
an  error  in  the  AFCRL  line  atlas.  In  Fig.  14(a).  the  relative  strengths 
appearing  in  the  AFCRL  line  atlas  for  the  strong  lines,  belonging  to  the  010 
to  030  band,  that  lie  between  2930  and  3000  cm"  have  been  represented  by 
vertical  lines.  There  is  a very  good  correlation  between  the  position  of 
these  lines  and  the  extra  features  in  the  computed  spectrum.  The  strength 
of  this  band  has  been  approximated  by  summing  the  strengths  of  all  lines 
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cm  to  3000  cm  . These  lines 
tures  in  the  calculated  spectrum. 
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in  the  atlas  for  this  band.  This  sum  has  the  value  1.46  X 10‘20  cm2  mole"1 

cm-1.  The  listed  strength  for  this  band  (Table  8,  Ref.  2)  is  7.99  X 10-23 

2 - 1 - 1 

cm  cm  mole  , but  J.  E.  A.  Selby  of  AFCRL  indicates  the  correct 

-23  2 -1  -1 

strength  is  1.  536  X 10  cm  mole  cm  . In  either  case,  the  lines  in  our 
copy  of  the  atlas  appear  to  be  too  3trong  by  about  3 orders  of  magnitude; 
this  is  confirmed  by  comparison  with  experimental  data  as  indicated  in  Fig. 
14(a).  Recent  recomputation  of  the  spectrum  with  a revised  atlas,  in  which 
the  intensity  of  this  band  was  corrected,  has  resulted  in  a spectrum  that 
agrees  much  better  with  the  experimental  data. 

The  band  strength  discrepancy  discussed  above  for  the  010  to  030  band 
is  by  far  the  largest  discrepancy  found  in  comparing  the  band  sums  from  the 
line  atlas  with  the  published  band  strengths.  In  view  of  the  vast  amount  of 
data  in  this  atlas,  it  is  understandable  how  such  an  error  could  occur.  Pres- 
ent experience  indicates,  however,  the  importance  of  comparing  computed 
spectra  with  experimental  spectra  at  sufficient  resolution  to  identify  indivi- 
dual spectral  features. 

D.  HOT  COz  PATH 

The  hot  CO  experimental  spectrum  indicated  by  the  solid  line  in  Fig. 

^ 10 

15  is  taken  from  the  work  of  F.  S.  Simmons.  This  spectrum  has  been 
replotted  by  hand  from  Fig.  6 of  Ref.  10.  The  calculated  spectrum  indicated 
by  the  dashed  line  was  obtained  by  conv  living  a triangular  instrument  function 
5 cm  1 wide  (FWHM)  with  the  high  rtrolution  spectrum  obtained  from  the 
AFCRL  line  atlas  with  the  use  of  the  INHOM  computer  program  for  the  exper- 
imental conditions  listed  in  Table  8.  Note  that  this  is  a radiance  spectrum 
rather  than  a transmittance  spectrum  such  as  those  presented  elsewhere  in 
this  section. 

The  excess  radiance  of  the  experimental  spectrum  compared  to  the 
calculated  curve  at  low  wave  numbers  is  due  to  transitions  not  included  in 
the  line  atlas  because  their  lower  levels  are  not  significantly  populated  under 
atmospheric  conditions.  This  defect  is  due  to  the  application  of  the  AFCRL 
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atlas  to  conditions  for  which  it  was  not  constructed.  The  excess  of  the 
calculated  radiance  over  the  experimental  at  high  wave  numbers,  on  the  other 
hand,  is  not  understood,  but  perhaps  represents  a problem  in  the  calibration 
of  the  experimental  setup.  Dips  in  the  radiance  in  the  calculated  spectrum, 
which  are  not  reproduced  in  the  experimental  spectrum,  may  be  due  to  the 
use  of  a spectral  slit  width  greater  than  5 cm"1  in  the  determination  of  the 
experimental  spectrum,  or  to  transitions  omitted  from  the  line  atlas. 

In  general,  in  this  spectral  region,  the  AFCRL  atlas  again  appears 
correct  for  CO^  except  for  the  intentional  omission  of  transitions  arising 
from  highly  excited  lower  levels.  These  omissions  can,  of  course,  lead  to 
serious  problems  in  the  wings  of  bands,  as  indicated  in  Fig.  15. 

E.  HOT  HzO  PATH 

The  experimental  spectrum  represented  by  the  solid  line  in  Fig.  16 
was  obtained  from  the  work  of  F.  S.  Simmons.  11  The  curve  illustrated 
v'as  obtained  from  a photograph  of  Fig.  47  in  Ref.  11  that  was  then  trans- 
formed into  a computer-compatible  format  by  the  method  described  in 
connection  with  the  high  resolution  atmospheric  CO^  path.  Here,  in  addition 
to  a shift  in  the  frequency  scale,  a 3-percent  compression  of  the  original 
wavelength  scale  was  required  to  allow  the  position  cf  the  major  features  of 
the  experimental  spectrum  to  agree  with  the  compu'ed  spectrum.  The  com- 
puted spectrum,  indicated  by  the  dashed  line,  was  calculated  by  convolving 
a 7.  5 cm"1  wide  (FWHM)  triangular  instrumental  function  with  the  high 
resolution  spectrum  obtained  with  the  INHOM  Program  from  the  AFCRL 
atlas  for  the  experimental  conditions  listed  in  Table  8. 

Agreement  of  the  experimental  and  computed  spectra  is  good  for 
frequencies  greater  than  about  3450  cm  ^ Around  4100  cm  there  is 
an  indication  of  the  effect  of  omitted  lines  from  highly  excited  states,  but 
the  difference  is  not  large.  Below  3450  cm"  the  comparison  is  very  poor. 


I 
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discussed  previously.  Several  of  these  lines  have  lower  energy  levels  of 
about  2000  cm-1,  so  that  increasing  the  temperature  from  300  K to  1000  K 
results  in  an  increase  of  about  3 orders  of  magnitude  in  the  Boltzmann  factor 
in  the  line  strength.  As  a result,  these  lines  are  strong  absorbers  in  this 
region.  Because  of  the  presence  of  these  incorrect  strong  lines,  no  conclu- 
sions can  be  drawn  about  the  applicability  of  our  version  of  the  AFCRL  line 
atlas  to  calculation  of  the  radiation  from  hot  l^O. 
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V.  SUMMARY 


• i program  and  associated  procedures  have  been  developed  to 
compute  high  resolution  spectra  for  inhomogeneous  optical 
paths,  including  the  conditions  present  in  a missile  plume. 

These  programs  and  procedures  are  outlined  in  this  report. 

• These  programs  and  procedures  have  been  used  with  the 
AFCRL  line  atlas  to  demonstrate  that  ignoring  the  correlation 
of  atmospheric  absorption  lines  with  source  emission  lines 

can  lead  to  errors  as  large  as  79  percent  in  the  computed  trans- 
mittance under  some  conditions. 

• The  AFCRL  line  atlas  needs  to  have  lines  added  to  it  from 
high.y  excited  states  to  provide  an  adequate  data  baBe  for  the 
calcalation  of  high  temperature  optical  properties  of  gases. 
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APPENDIX 

USER'S  GUIDE  TO  IN  10M 


1.  INTRODUCTION 

The  understanding  of  mary  atmospheric  optical  phenomena  can  be 
studied  by  modeling  these  optical  effects  with  computer  calculations.  The 
INHOM  Program,  which  is  described  here,  may  be  used  to  calculate  the 
transmittance  and  radiance  of  an  optical  path  through  an  atmosphere  which 
may  or  may  not  be  homogeneous.  This  program  is  structured  modularly 
so  that,  by  replacing  various  modules  in  the  program,  almost  any  physical 
process  that  can  be  incorporated  into  a radiative  transfer  problem  can  be 
acc  ommodated . 

The  programs  have  been  written  principally  in  FORTRAN  — IV  com- 
patible with  the  Control  Data  Corporation  (CDC)  FTN  compiler  in  use  on  the 
CDC  7600  computer  at  The  Aerospace  Corporation.  A few  routines,  which 
must  be  very  efficient  or  must  perform  functions  not  readily  programmable 
in  FORTRAN,  have  been  written  in  the  CDC  assembly  language  COMPASS. 

The  first  application  of  this  program  has  been  to  compute  infrared 
atmospheric  transmittance  and  radiance  under  conditions  typical  of  the  lower 
atmosphere  where  the  assumption  of  local  thermodynamic  equilibrium  is  valid 
and  the  individual  spectral  lines  are  pressure  broadened.  The  basic  spectro- 
scopic input  to  these  calculations  is  the  atlas  of  line  parameters  prepared 

under  the  coordination  of  R.  A.  McClatchey  at  the  Air  Force  Cambridge 

2 

Research  Laboratories,  The  format  of  this  input  data  is  described  in  para- 
graph 3.  The  output  of  the  program  is  a printed  listing  of  computed  trans- 
mittance and  radiance  together  with  a listing  of  various  intermediate  results. 
Trarsmittance/radianco  results  are  also  available  in  computer  compatible 
format  and  may  be  saved  on  magnetic  tape.  The  format  of  this  output, 
described  in  paragraph  4,  is  compatible  with  the  input  requirements  for  the 
spectral  analysis  program  SPECPST  already  available  (see  Ref.  3,  Section  V). 

I 

I 


Preceding  page  blank 
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Details  of  the  specific  calculation  to  be  performed  are  controlled  by  computation 
step  cards  described  in  paragraph  2.  A sample  run  is  described  in  paragraph  5. 

2.  INHOM  PROGRAM  CONTROL 

a.  INTRODUCTION 

The  computations  performed  by  the  INHOM  Program  are  controlled  by  a 
series  of  computation  step  cards.  These  cards  are  read  as  input  data  by  the 
computer  program  after  it  has  been  called  into  execution. 

Each  computation  of  a spectrum  is  initiated  by  a SPECTRUM  step  card 
and  results  in  the  computation  of  a single  spectrum.  The  physical  parameters 
for  which  the  spectrum  is  computed  are  specified  by  parameters  on  the 
SPECTRUM  step  card  and  the  information  available  from  other  preceding 
step  definition  cards.  All  information  input  by  step  cards  other  than  a 
SPECTRUM  step  card  is  retained  unchanged  by  the  program  until  it  is  super - 
ceded  by  the  occurrence  in  the  input  stream  of  the  same  step  definition  with 
different  parameters,  in  this  way,  it  is  possible  to  calculate  spectra  under 
varying  physical  conditions  by  respecifxung  only  those  variables  that  change 
rather  than  all  the  variables  influencing  the  calculation. 

b.  STEP  DEFINITION  FORMAT 

All  computation  step  cards,  except  as  explicitly  noted  in  the  detailed 
step  descriptions,  have  the  following  properties  in  common. 

(1)  The  cards  are  free  form  in  the  sense  that  there  is  no  rigid 
requirement  that  specifications  begin  in  a particular  card 
column  or  that  they  have  a specific  fixed  fiel  i length.  The 
beginning  and  length  of  fields  are  specified  explicitly  on  the 
input  cards  by  certain  symbols  and  conventions. 

(2)  The  f.pecification  of  each  step  must  begin  on  a new  card. 

In  some  cases,  this  may  extend  for  more  than  one  card. 

(3)  Spelling  is  important.  Only  the  combin?tions  of  letters 
and  numbers  listed  in  the  detailed  step  descriptions  are 
recognized.  All  other  combinations  will  result  in  error 
messages  to  the  user  from  the  program.  No  blanks  may 
be  embedded  wit  un  name!'  and  numerical  values  since  a 
blank  is  used  as  a field  delimiter. 
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(4)  The  first  word  on  the  card,  called  a step  definition  word, 
indicates  to  the  program  what  is  to  be  done. 

(5)  The  order  in  which  steps  are  specified  preceding  a SPECTRUM 
step  is  usually  not  important.  When  a SPECTRUM  step  is 
encountered,  it  computes  a spectrum  based  on  all  the  information 
specified  to  that  point  in  the  input  stream. 

(6)  For  some  steps,  numerical  parameters  are  required.  These 
are  input  by  a routine  similar  to  an  internal  FORTRAN 
NAMELIST  routine. 

(a)  The  beginning  of  the  parameter  list  is  indicated  by  a 
left  parenthesis 

(b)  The  parameters  are  specified  in  the  format  name  = 
value.  In  the  step  descriptions  when  the  value  is 
indicated  by  nnn,  a fixed  point  number  is  required. 

When  the  value  is  indicated  by  fff,  a floating  point 
number  is  required.  If  the  decimal  point  is  omitted 
from  a floating  point  number,  it  will  be  assumed 
immediately  to  the  right  of  the  last  digit. 

(c)  Specifications  are  separated  by  commas  and  terminated 
with  a slash,  /,  or  right  parenthesis,  ). 

(d)  No  blanks  may  be  embedded  in  names  or  values,  but 
they  are  permitted  elsewhere. 

(e)  Parameter  definition  must  begin  on  the  first  card  of  the 
step  or  definition  and,  unless  stated  otherwise,  must 
be  completed  on  that  card. 

(f)  The  parameters  may  be  specified  in  any  order, 
c.  DETAILED  COMPUTATION  STEP  DEFINITION 

TITLE 

The  contents  of  Columns  11  through  80  of  his  card  are  concatenated 
with  the  date  on  which  the  run  was  made.  The  resulting  80 -character - 
long  string  is  used  whenever  a title  is  needed  in  the  printed  listing. 

This  string  is  also  used  to  label  the  output  spectrum  put  onto  magnetic 
tape.  A title  remains  unchanged  until  another  TITLE  step  is  encountered. 

STOP 

This  step  ends  the  run.  Any  step  control  cards  remaining  in  that 
section  of  the  input  file  are  ignored. 
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atmosphere 

This  step  specifies  the  input  atmo'phe're  th  "“iL^tmos^he"^'^ 
K'ipecmed  i^rS^  o/ temperature  ^ pressure. 

SSSW^  i;=r^to  «nf  - 

sfi'u=SHSiH5=i:. 

The  atmosphere  specified  by  this  step  rema 
next  ATMOSPHERE  step  is  encountered. 


Step  Card  For  aat: 

ATMOSPHERE 


Inone  ) none 

direct}  direct 

PRINT  \ PRINT 

The  atmosphere  is  specified  on  the  ca r d s d'on ^the^i ATMOSPHE RE 

card  in  a format  depending  on  the  option  !''“d  d,h"  direct  option  will 

word  PRINT  appears,  the  input  atmosphere  v,  x * P 
put  listing;  otherwise,  it  will  not  be  printe  . 

(1)  Altitude  Profile  Option 

The  first  card  following  the  ATMOSPHERE  card  contains  in 
Columns  and  6 the  number  of  layers  in  the  atmosphere  This 
°r0  the  number  of  cards  to  follow.  A maximum  of  81  layers 

may  he  specified  The ^ ,*« ^^ci^g 

:ri;ndtiLr,trmo.;Ler“c\evIus  indicated  in  the  foil lowing 
listing.  All  variables  are  floating-point  real  number  . 

Altitude  Profile  Atmosphere  Input  Level  Specification  Format 


Card  Col  imna 


Physical  Quantity 


Units 


1 

- 6 

Altitude 

kilometers 

millibars 

7 

- 16 

Pressure 

17 

- 22 

Temperature 

Kelvin 
. 3 

23 

- 32 

Water  Vapor  Density 

gm / m 

3 

33 

- 42 

Ozone  Density 

gm/ m 

43 

- 56 

Aerosol  Number  Density 

particle  s/ cm 
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(2)  Direct  Input  Option 

The  cards  following  the  ATMOSPHERE  card  contain  the  path 
length,  molecular  concentration,  temperature,  and  pressure 
in  each  segment  of  the  atmosphere  through  which  the  optical 
path  passes.  The  input  for  each  segment  is  a free  form 
parameter  list.  Each  segment  specification  begins  on  a new 
card  and  is  indicated  by  a left  parenthesis,  (.  Specification 
of  that  level  continues  until  a right  parenthesis,  ),  is  encountered. 
Specification  of  a single  segment  may  be  continued  over  several 
cards  by  starting  a new  card  anywhere  a comma  is  allowed  in 
the  free  form  parameter  list.  Segments  are  assumed  to  be 
specified  in  the  order  in  which  they  will  occur  along  the  optical 
path.  The  first  left  parenthesis  introduces  the  first  segment 
of  the  optical  path.  This  segment  specification  is  terminated 
by  the  first  right  parenthesis.  The  next  segment  begins  on  the 
next  card  with  a left  parenthesis,  and  so  on.  The  direct  input 
option  is  terminated  by:  (a)  completing  a segment  specification 
with  a right  parenthesis;  and  then  (b)  following  this  with  one  card, 
which  may  have  any  chararte~s  on  it  except  a left  parenthesis. 

The  variable  names  recognized  and  the  units  associated  with 
them  are  as  follows: 


Variable  Name 


F y sical  Quantity 


Units 


p 

= 

fff 

Pressure 

Millibars 

T 

= 

fff 

Temperature 

Kelvin 

WH20 

= 

fff 

Water  Concentration 

molecules/cm 

WC02 

r 

fff 

CO^  Concentration 

molecules/cm 

W03 

= 

fff 

O^  Concentration 

molecules/cm 

WN20 

= 

fff 

N^O  Concentration 

molecules/cm 

WCH4 

r 

fff 

CH„  Concentration 
4 

molecules/cm 

WCO 

= 

fff 

CO  Concentration 

molecules /cm 

W02 

= 

fff 

Oz  Concentration 

molecules/ cm 

WHAZE 

= 

fff 

Aerosol  Concentration 

2 

particles/cm 

DX 

= 

fff 

Path  Length  in  Segment 

kilometers 

COMPONENTS 


This  controls  which  gases  are  assumed  to  exist  in  the  atmosphere. 
If  gases  are  turned  off  here,  the  lines  are  never  removed  from  the 
atlas  and  less  computing  time  is  required  than  if  the  concentration 
is  set  to  zero. 
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Format: 


COMPONENTS  name  j^J^j 

where  name  is  one  of  fVr*  following  words: 

H20 
C02 
03 
N20 
CO 
CH4 
02 

HAZE 
ALL 

The  OFF  following  one  of  these  words  tells  the  program  to  omit 
that  component  from  the  computation. 

ON  following  one  of  these  words  includes  that  component  in  the 
computation. 

ALL  turns  all  of  the  components  off  or  on  depending  on  the  word 
following. 

Words  must  be  separated  by  spaces,  commas,  or  left  parentheses,  (. 
Initially  all  components  are  on. 

GEOMETRY 

This  specifies  the  optical  path  followed  by  the  radiation  through  the 
atmosphere.  This  step  will  convert  an  atmosphere,  specified  by 
the  altitude  profile  option,  into  the  parameters  required  for  calcula- 
tions along  the  optical  path.  If  the  ATMOSPHERE  DIRECT  option 
is  used,  there  is  no  conversion  to  be  performed  and,  therefore,  data 
input  with  a GEOMETRY  step  are  ignored.  The  actual  conversion 
takes  place  when  the  computation  is  initiated  by  the  SPECTRUM  step, 
so  the  order  of  ATMOSPHERE  and  GEOMETRY  steps  in  the  input 
stream  is  unimportant.  The  GEOMF  TRY  specified  will  continue  to 
be  used  until  another  GEOMETRY  card  is  encountered.  Unless  the 
ATMOSPHERE  DIRECT  step  is  specified,  at  least  one  GEOMETRY 
step  must  precede  the  first  SPECTRUM  step. 

Refraction  is  ignored.  The  earth  is  assumed  to  be  spherical  with  a 
radius  of  6378.  16  km.  Impossible  geometric  paths  are  flagged.  All 
possible  gt  ometric  paths,  including  zenith  angles  greater  than  90  deg, 
are  believe!  to  be  handled  correctly. 


(OFF 
| ON 
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Format; 


GEOMETRY  (Parameters) 


Parameters : 


BEGALT 

FINALT 

ZA 

NLEVEL 

LNGPTH 


= fff  is  the  beginning  altitude  in  km  for  the  path. 

= fff  is  the  final  altitude  in  km  for  the  path. 

fff  is  the  zenith  angle  of  the  optical  path  at 
the  beginning  point  in  degrees. 


nnn  is  the  number  of  segments  into  which  the 
path  is  broken.  A maximum  of  81  is  allowed. 


i indicates  for  ZA  > 90  deg  that  the  longer  of 
the  two  possible  paths  is  to  be  considered.  Any 
other  value  results  in  considering  the  shorter 

f°T  ZA  “9°  deg’  this  variat>le  is  ignored 
I he  default  value  is  zero. 


PRINT 


(on  ) 

/ 

(anything  else! 

type  J 

This  controls  the  printing  of  auxiliary  information.  The  print 
options  remain  in  force  until  overridden  by  another  PRINT  step. 

Format: 

l _ \ 

ON  ) 

anything  elsej 

ni?6  inKi1CaltS  the  information  desired  for  auxiliary  output  and 

i're;ir or icular  prin,in8  op,i°"-  -k  ««• 

The  types  recognized  are: 

ADDLINES  Every  time  additional  lines  are  extracted 
from  the  line  alias,  they  are  listed. 

W MATRIX  Before  the  spectrum  computation  is  initiated, 
the  matrix  of  environmental  parameters 
specifying  the  segments  along  the  optical  path 
will  be  printed. 

SHAPE  The  line  strength  and  line -broadening  parametei 
to  be  used  for  each  segment  are  printed. 

DIVY  The  altitudes  of  the  ends  of  the  segments  used  fc 

converting  an  atmosphere  to  a segmented  path  a 
printed. 
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These  print  options  are  initially  set  to  OFF  (no  printing). 

It  should  be  noted  that  some  other  auxiliary  printing  options  are 
controlled  by  parameters  in  certain  step  definitions. 

SPECTRUM 

This  step  initiates  the  computation  of  a spectrum.  It  is  the  one  step 
that  is  required  for  every  spectrum. 

Format: 

SPECTRUM  (Parameters) 

Parameters : 


The  first  four  parameters  must  be  jpecified  on  every 
SPECTRUM  step  card.  The  other  parameters  are 
optional  and  will  be  reset  to  their  default  v-.lues  for 
every  computation  unles s otherwise  specified. 


BEG 

END 

DNU 

BOUND 


fff  is  the  beginning  frequency  for  which  the 
spectrum  is  desired  (cm‘^). 

fff  is  the  ending  frequency  for  which  the 
spectrum  is  desired  (cm'*). 

fff  is  the  step  size  by  which  frequency  is 
stepped  from  BEG  to  END  (cm-1). 

fff  is  the  maximum  distance  away  from  a 

frequency  for  which  lines  are  considered 
to  contribute  (cm'*), 


IPNT  = 1 

ITRN  * 0 


causes  the  radiance  and  transmittance  to 
be  printed  out.  Other  values  cause  no  listing. 
The  default  is  no  listing. 

causes  only  the  transmittance  spectrum  to  be 
computed.  This  saves  time  and  makes  each 
output  spectrum  record  contain  twice  the 
frequency  interval.  The  default  is  zero  and 
so  both  the  radiance  and  transmittance  are 
computed. 


If  both  radiance  and  transmittance  are  com- 
puted, the  spectrum  type  (Table  A-3,  word  9) 
is  LINE  SPECT.  If  only  transmittances  are 
computed,  the  spectrum  type  is  SPECTRUM. 
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3. 


LINE  ATLAS  FORMAT 


The  spectroscopic  data  required  for  the  calculation  of  transmittance  is 
stored  on  magnetic  tape  in  the  form  of  80-character  punched-card  images. 
Each  card  image  corresponds  to  one  line,  and  up  to  40  lines  ~ard  images) 
are  stored  in  each  block  on  the  tape.  The  first  word  of  each  block  contains 
an  integer  indicating  the  number  of  lines  stored  in  that  block.  The  lines  are 
organized  in  an  ascending  order  of  frequency.  For  convenience  in  gaining 
iccess  to  the  data,  different  spectral  regions  are  separated  into  different 

files.  These  spectral  regions  and  their  associated  file  numbers  are  shown  in 
Table  A - 1 . 


Table  A-l.  Frequency  Range  Covered  by  Each 
File  on  the  Line  Atlas  Tape 


File 

Frequency  Interval 
(cm-1) 

1 

0 - 500 

2 

500  - 1000 

3 

1000  - 2000 

4 

2000  - 5000 

5 

5000  - 10,000 

Spectroscopic  parameters  needed  for  the  calculation  are  listed  in 
Table  A-2  along  with  the  fields  on  the  card  employed  and  the  FORTRAN 
format  associated  with  the  field.  For  further  details,  consult  Ref.  2. 

At  present,  about  110,000  lines  are  available  in  the  atlas.  These 
include  seven  gases  listed  in  Table  A-2. 


Table  A-2.  Parameters  in  Spectroscopic  Line  Atlas 


Characters 


Parameter 


Line  center  frequency  cm 

Line  strength  at  296  K cm"  1 /(molecules  cm  2) 

Line  width  parameter  cm  * /atmosphere 
at  296  K 

.1 

Energy  level  of  lower 
state  in  transition 

Quantum  numbers 

Date 

a 

Isotope 

Molecule  indicated 
as  follows: 


Molecule  Code 


The  detailed  code  depends  on  which  molecule  is  involved  (see  Ref.  2) 


COMPUTER  COMPATIBLE  OUTPUT  FORMAT 


■ 


! 

I 

i 

* 

§ 


5Br 


The  computed  spectra  are  available  in  computer  compatible  format  on 

file  TAPE30.  The  format  of  this  file  is  compatible  with  the  input  to  SPECPST, 

which  is  the  spectral  analysis  program  that  is  part  of  the  Fourier  spectroscopy 

3 

computer  package  available  at  Aerospace. 

When  the  file  is  put  on  tape,  it  is  7 -track  800 -bpi  odd  parity.  A 
spectrum  may  contain  one  or  more  logical  records.  The  format  of  each 
logical  record  is  shown  in  Table  A-3.  The  length  of  every  logical  record 
in  a single  spectrum  is  the  same,  but  the  logical  records  in  different  files 
on  the  same  tape  may  be  of  different  lengths. 

5.  SAMPLE  PROBLEM 

A sample  computer  run  is  shown  in  Figs.  A- 1 and  A-2.  In  Fig.  A-l, 
a sample  input  deck  is  illustrated  that  uses  the  altitude  profile  method  of 
specifying  the  atmosphere  through  which  the  radiation  passes.  Figure  A-2 
shows  the  computed  results  produced  by  INHOM  on  the  basis  of  the  input 
deck  illustrated  in  Fig.  A-l. 


a 


I 
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Table  A -3.  Logical  Record  Format  l'or  Spectrum  Output  Tapes 


% ord 
No. 

Type 

Description 

1 

Integer 

File  number.  This  begins  with  1 for  the  first  file  on 
the  tape. 

2 

Integer 

The  number  of  logical  records  in  this  file. 

3 

Integer 

The  logical  record  number,  K,  of  this  record.  The  first 
logical  record  is  No.  1. 

4 

Floating 

The  frequency  increment  between  data  values  in  this  file. 

5 

Floating 

The  frequency  corresponding  to  the  firs  ue  in  the 

spectrum. 

6 

Floating 

The  maximum  frequency  for  which  a value  is  specified  in 
the  spectrum. 

7 

Floating 

The  minimum  value  in  the  spectrum.  (For  LINE  SPECT 
spectra,  this  is  zero.  ) 

8 

Floating 

The  maximum  value  in  the  spectrum. 

For  LINE  SPECT  spectrum,  it  is  the  maximum  value  in 
the  radiance  spectrum. 

9 

Display 

Spectrum  type. 

SPECTRUM  indicates  a transmittance -only  spectrum 
obtained  by  setting  ITRN  * 0 in  the  SPECTRUM  step. 

LINE  SPECT  indicates  the  output  consists  of  a 
transmittance  and  radiance  spectrum. 

10-17 

Display 

The  title  used  to  label  the  printout  when  the  tape  was 
originally  generated. 

18 

Integer 

NWRDS,  the  number  of  remaining  words  in  this  logical 
record,  all  of  which  are  data  words. 

19  to 

NWRDS+1 8 

Floating 

The  values  of  the  spectrum.  Details  depend  on  the 
spectrum  type.  On  multirecord  files,  it  is  the  values 
(K-1)*N  + t to  K*N  where  N depends  on  the  spectrum 
type  (defined  below)  and  K comes  from  Word  3. 

For  SPECTRUM  type  spectra,  N = NWRDS.  All  of  the 
values  are  spectral  values  in  ascending  order  of 
frequency. 

For  LINE  SPECT  type  spectra,  N = NWRDS/2.  Words 
1 to  N of  the  values  section  are  the  transmittance  values 
in  ascending  order  of  frequency.  Words  N + l to  NWRDS 
are  the  radiance  values  in  ascending  order  of  frequency. 

i«313F*C3  300.0 
1*0  9.040t*02  294.0 

|*S  S*?|2f*°2  288*0 
3.0  7. 150E+02  2 H.C 
f»*0  6.330E+02  2’7.0 
J*2  ? *1 29  0 2 270.0 

7*-  £*li?f+92  264.0 
7*  4 . 320  E*r  2 267 . C 

8*0  3.M0E702  250.0 

.?•$  5 *5 29|^°  2 244.0 
1 • 2.96JE*  2 2 37  , C 

i*e2E*02  230.0 
*2*0  2.130E+02  224. C 
13. t 1.920^4.02  217.0 

14.0  1.560E+02  210.0 

}£•£  f *?29f + 0 2 204.0 
16*0  1 .1 10E+02  197.0 

}!•  3*H?E+C1  115.c 

18.  7 . 3 9 0 1 19  9.' 

?n*n  203. C 

21.  4.9oCE>i.  1 211,1 

I?*n  7‘cnJlf*01  215.0 

23.0  3.500E701  217.0 

24.0  3.000E>01  219.0 
2^.C  2 .570E4-C  1 221.0 
H*9  i*220Ef01  232. C 
.35*u  S*000t*00  2-3,0 

!rH*n  ?*-I2E*00  2*4. C 
crt*2  i*ig0E4-0C  265.  t 
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Fig.  A-?.  Spectrum  Computer!  and  Listing  Gene  rated  by  the  INHOM  Program 
in  Response  to  the  Input  Deck  Shown  in  Fig.  A-i  (Continued) 


uSSSSS=;SSS§SSS|SSS§SSSS2^SS 

^WhiWWItttiMyiiMiii.iiiii.ii.M  a *,  * * L*  ♦ ♦ ♦ 


g^s^assr'SssSsS 


3hhhhhh< 


p*BBBB$BBBBBBB** 


fY 

UJ 


STSSSS^  I 

s-S^^rs2;c2-s5^£s;^^sss^5;ss  § 


U^gOOQOOO  JQOOOOOQ 

g^5S5SSS23SS5SSsS^i2^S5 


'OOOOiOq  JOO 
— OOOOOOjflJOO 
444444.*4444  , 4 1444 


«\J 

u 


aR8s»ss=SsS55^^^ 



gSSiS^ 

rSSSSSSS2-?J.c5  Q 


'i  ill'  ••••••  .*.  • • • • • • . vrrrrrrrrlu^ 


-•(/» 

„ r 

^OO^HOOO  7 

ooroooo  o 

c 

UJUJ 
-it'' 

- o ir  ^^«f^rrr  •■‘(r-rk  tft^'ir  'f  x » cr a-  rr»-« 


*444*444444444  4444  444  7*  + 1**1 

?rlJw%^!i."^i"Ll'IUJU  tiiU-UJl'JU^  ItJUJUj  UJ  UJ  i . 
iZ-LiT^T  <-t~K  •..fv»'rf>(-i».  o^_.ir>oifvr. 


L Lrir*.rtn.rZ”’t.  ^ v 

tlinif  inininininf!^^^^^^^^^x^’'fl'*’^^’2^l25tnZ. 


- -u  - 


£ 

<13 

(1 

60 

0 „ 
Ul  OJ 

,_,  T3 
2 £ 
O y 
X C 

XU 

<D 

£7 

-O  . 
•O.M 
-2  U, 

2 C 
a»  ■H 
C C 

O o 

?“ 

£ X 

m o 

• H 4/ 

JQ 

' V *J 

g & 

<4  g- 
T?  ‘-, 
^ <u 

3 -c 
a" 
£2 


o 

U 


o 
m 

£ 0 
p a 

>-»  cn 
-*-»  o 

o*  c 

w .1: 


<\) 

I 

< 

60 
• H 

u, 


REFERENCES 


S.  J.  Young,  Band  Model  Calculations  of  Atmospheric  Transmittance 
for  Hot  Gas  Line  Emission  Sources,  TR- 0075(5647)- 1 . The  Aerospace 
Corporation,  El  Segundo,  California  (in  preparation,  June  1974). 

R.  A.  McClatchey,  W.  S.  Benedict,  S.  A.  Clough,  D.  E.  Burch, 

R.  F.  Calfee,  K.  Fox,  L.  S.  Rothman,  and  J.  S.  Garing,  AFCRL 
Atmospheric  Absorption  Line  Parameters  Compilation,  Environ- 
mental Research  Paper  434,  AFCRL- 'iR-73-0096,  "Air  Force 
Cambridge  Research  Laboratories,  Mass.  (26  January  1973). 

C.  M.  Randall,  User's  Guide  to  Computer  Programs  for  Fourier 
Spectroscop/,  TR-0073(9260-01 )- 8,  The  Aerospace  Corporation, 

El  Segundo,  California  (30  November  1972). 

G.  Herzberg,  Molecular  Spectra  and  Molecular  Structure,  Van 
Nostrand  Reinhold  Company,  New  York,  New  York  (1945),  p.  501. 

C.  B.  Ludwig,  W.  Malkmus,  J.  E.  Reardon,  and  J.  A.  L.  Thomson, 
Handbook  of  Infrared  Radiation  from  Combustion  Gases,  NASA 
Special  Publication  NASA  SP-3080  (1973). 

R.  A.  McClatchey,  R.  W.  Fenn,  J.  E.  A.  Selby,  F.  E.  Volz,  and 
J.  S.  Garing,  Optical  Properties  of  the  Atmosphere  (Third  Ed.), 
Environmental  Research  Paper  No.  411,  AFCRL- 72-0497 , Air  Force 
Cambridge  Research  Laboratories,  Mass.  (24  August  1972). 

J.  E.  A.  Selby  and  R.  A.  McClatchey,  Atmospheric  Transmittance 
from  0.25  to  28.5  pm:  Computer  Code  LOWTRAN2,  Environmental 
Research  Paper  No.  427,  AFCRL-72-0745,  Air  Force  Cambridge 
Research  Laboratories,  Mass.  (29  December  1972). 

D.  E.  Burch,  D.  A.  Gryvnak,  and  R.  R.  Patty,  Absorption  by 
CO2  Between  3100  and  4100  cm-  1 (2.44  - 3.22  Microns),  U-4132, 

Aeronutronic  Division  of  Philco-Ford,  Newport  Beach,  California 
(30  April  1968). 

D.  E.  Burch,  D.  A.  Gryvnak,  and  R.  R.  Patty,  Absorption  by  H?Q 
Between  2800  and  4500  cm*'  (2.  7 Micron  Region),  U-3202, 
Aeronutronic  Division  of  Philco-Ford,  Newport  Beach,  California 
(30  September  1965). 


Preceding  page  blank 


-83- 


10. 


\ 


F.  S.  Simmons , H.  Y.  Yamada,  and  C.  B.  Arnold,  Measurements 
of  Temperature  Profiles  in  Hot  Gases  by  Emission- Absorption 
Spectioscopy,  University  of  Michiean  Report  WRL  8962-lft-Fr 
NASA  CR-7?491,  Willow  Run  Laboiatories  (April  1969). 

11.  F.  S.  Simmons,  C.  B.  Arnold,  and  D.  H.  Smith,  Studies  of  Infrared 
Radiative  Transfer  in  Hot  Gases  1:  Spectral  Absorptance  Measurc- 
ments  in  the  2.  7|ji  H2O  Bands,  University  of  Michigan  Report 

4613-91-T,  Willow  Run  Laboratories  (August  1965). 

- 

. 

. 

1 


.1 


j 

j 


1 


f 


-84- 


* 


' 


laboratory  operations 

Th.  Laboratory  Operation,  of  The  Aero.p.ce  Corporation  i.  conducting 
experimental  and  theoretical  inve.tigation.  n.cee.ary  for  th.  evaluation  and 
application  of  acienHflc  advance,  to  new  military  concep;.  and  eyeteme.  Ver- 
eatiiity  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
pereonnel  in  dealing  with  the  many  problem,  encountered  in  the  nation-.  rapidl> 
developing  .pace  and  mia.ile  .y.t.m.,  Experti.e  in  the  late.t  .cientific  devel- 
opment. i.  vital  to  th.  accomplishment  of  task,  related  to  the.e  problem..  The 
laboratorle.  that  contribute  to  thi.  research  are: 


, Aerophy.ic.  Laboratory:  Launch  and  reentry  aerodynamic,  heat  tram 
f!r’  r®*nt*V  phy.ic.,  chemical  kinetics,  .tructural  mechulci.  night  dynamic. 
atmo.ph.ric  pollution,  and  high-power  ga.  la.er.,  dynamic, 

„l— l.^7.-e/inl-*tryL*n<l  P|Y*1C*  L>bor>torr:  Atmo.pherlc  reaction.  and  atmo.- 

of  e x'c ited*  « p e c io  7n  r cick f ' V °" * Pv"Ut*d  •tmo*P*'«>-e.,  chemical  reaction. 

, excited  specie.  in  rocket  plume.,  chemical  thermodynamic  nla.m.  anH 

ind'c^H  r*4C.tl0n,‘  ‘b.mi.try,  propulsion  cne^.trV,  <p.c  uum 

cn.mv.  mat  . , °d  lubrication  and  .urf.ce  phenomena,  photo- 

maging,  atmo.pheHc  pollution;  millimeter  wave  and  far-infrared  technology. 

Material.  Science.  Laboratory:  Development  of  new  material.;  metal 

matrix  composite,  and  new  form,  of  carbon;  test  and  evaluation  of  xrauhite 
and  ceramic,  in  reentry;  spacecraft  material,  and  electronic  -omponent.  in 

rcionr.nH  fP«n<  en1Wr°n™rnti  »PPli cation  of  fracture  mechanic,  to  .tree  cor- 
rovion  and  fatigue-induced  fracture*  in  structural  metal*. 

Hon  Phy,ic*  Laboratory:  Atmo.pheric  and  ionospheric  phy.ic.,  radia- 

atmosphere,  d.n.ity  and  compo.ition  of  the»tmo.£h^  iuror.e 
and.irglow,  magneto.pheric  phy.ic,  co.mic  ray.,  generation  and  p^opaction 
of  pla.ma  waves  in  the  magneto.phere;  solar  phy.ic.,  studio  of  solar  m/ineHc 

mVeneHc  *1*  *,trono"'y-  x-r»Y  astronomy;  the  effect,  of  nuclear  explo.ion. 
magnetic  .torm.,  and  .olar  activity  on  the  earth1,  atmo.phere  iono.ohere  and 
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